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ABSTRACT 

We present extensive optical spectroscopy of the early-type magnetic star HD 191612 
(06.5f?pe-08fp). The Balmer and He i lines show strongly variable emission which is 
highly reproducible on a well-determined 538-d period. Metal lines and He 11 absorp- 
tions (including many selective emission lines but excluding He 11 A4686A emission) are 
essentially constant in line strength, but are variable in velocity, establishing a double- 
lined binary orbit with Poib = 1542d, e = 0.45. We conduct a model-atmosphere anal- 
ysis of the primary, and find that the system is consistent with a ~08 giant with a 
~B1 main-sequence secondary. Since the periodic 538-d changes are unrelated to or- 
bital motion, rotational modulation of a magnetically constrained plasma is strongly 
favoured as the most likely underlying 'clock'. An upper limit on the equatorial rota- 
tion is consistent with this hypothesis, but is too weak to provide a strong constraint. 



INTRODUCTION 



Peculiarities in the spectrum of the early-type star 



HD 191612 were first noted by Walborn (19731; it remains 
one of only three known Galactic examples of the Of?p 
class[^ this designation indicating C ill A4650 in emission 
with comparable strength to N ill A4640. Renewed inter- 
est followed the discovery of recurrent spectral variabil- 
ity between spectral types 06-07 and 08 (with correlated 



1 The others are HD 108 & HD 148937; the S MC Of?p stars 
AzV 220 & 2dFS 936 were found subsequentl y | |Walborn et al.| 



changes in the unusual emission-line features; [Walborn et al.| 
20031, which Walborn et al. (20041 showed to be consistent 
with a ~540-d period identified in Hipparcos photometry 
( Koen fc Eyer|2002|[Naze|2004| ). 



From timescale arguments, [Walborn et al.| ( |2004[ ) sug- 
gested that the 'clock' underlying the variability was most 
probably a binary orbit. New light was cast on this issue by 
[Donati et aL] (pb06a), who found HD 191612 to be only the 
second 0-type star known to possess a magnetic fieldrl Al- 
though their observations sampled only a single epoch, they 
were none the less able to estimate a polar field strength of 



2000[|Massey fc Duffy|2001[ [Evans et al.|2004 l 



2 The first was 6^ Ori C; Donati et al. I 2002 I 
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~1.5kG (from an observed line-of-sight field of ~220G, by 
assuming a dipole field), and made the case for 538-d rota- 
tional modulation, arguing that the field itself could easily 
be responsible for the implied slow rotation (through mag- 
netic braking). 

Of observational necessity, the magnetic, rotational- 
modulation model is as yet poorly constrained, and the dis- 
cussion of XMM- Newton spectroscopy by [Naze et al.| ( [2007[ ) 
emphasizes a number of discrepancies with the X-ray be- 
haviour expected in the simplest version of this scenario. 
On the other hand, the alternative orbital model has not 
been subject to any strong tests (arguably, even a strict, co- 
herent spectroscopic periodicity has yet to be demonstrated 
robustly). Here we present the results of an extensive cam- 
paign of optical spectroscopy, carried out in an attempt to 
shed light on these issues. 



2 OBSERVATIONS 

The major part of our campaign was conducted during the 
2004 and 2005 observing seasons. Because of the ~18-month 
variability timescale, scheduled observations at common- 
user facilities were generally impractical, and our observa- 
tions were obtained through service programmes; by taking 
advantage of telescope time awarded to other scheduled pro- 
grammes; and by exploiting the goodwill of colleagues. The 



Table 1. Ha variability: best-fit parameters for the arbitrary 
functional form described by eqtn^ 



main dataset, summarized in Table Al (166 digital observa- 
tions spanning 17 years), is therefore quite heterogeneous. 
None the less, the spectra can be conveniently character- 
ized by wavelength range ('red', including Ha 6563A; 'blue', 
generally including at least the ~4400-4700A region) and by 
resolution ('high', -R>4 x 10*; 'intermediate', R>4 x 10^; 
and 'low'). With a few exceptions, the spectra are generally 
reasonably well exposed, with signahnoise ratios typically 
approaching ~100. 

All spectra have been put on a heliocentric velocity 
scale, and all the Ha spectra have been corrected for telluric 
absorption by division, in topocentric space, with an appro- 
priately scaled and smoothed telluric 'map' constructed from 
high-quality, high-dispersion echelle spectra. (Because the 
telluric lines are unresolved, direct scaling in optical depth 
is impossible; we scaled in observed intensity, but checked 
that scaling in observed, pseudo-optical depth gives negligi- 
bly different results.) 



3 THE 538-D PERIOD: Ha VARIABILITY 



As already noted by |Naze et al.|p007[ ), the spectral lines can, 
for the most part, be separated into two groups according 
to variability characteristics: the absorption lines of metals 
and of He ll show, at most, small changes in line strength, 
while the hydrogen and He I lines show large equivalent- 
width changes. Ha shows the largest-amplitude variations 
of any spectral feature followed in our campaign, by a com- 
fortable margin. This, and the rather extensive temporal 
coverage of the red-region spectra, mean that it is the most 
useful feature for investigating the periodicity of the spec- 
troscopic variability. 
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Even casual inspection shows systematic variations in Ha 
that repeat on a period close to the ~540-d timescale found 
in Hipparcos photometry (over only about 2 cycles) . In order 
to quantify that period and its uncertainty, we fit an ad hoc 
analytical function to the equivalent-width measurements; 
we find that a truncated gaussian. 



Wx{(t>o,) = Wo-Aexp 



2al 

4> 



-00 < (j>a < -f 00 



00 < < 1 



(1) 

gives a good match to the observations. Although this func- 
tional form is arbitrary, it does provide a useful characteriza- 
tion of the data; the best-fit parameters are listed in Tablejl] 
and lead to the ephemeris 

0Q = (t - JD 2453415.2)/537.6d (2) 

where phase zero corresponds to peak Ha emission]^ 

The phase-zero epoch is chosen to be close to the median 
date of the observations (to minimize the formal error onto), 
and so, because the density of observations has increased 
with time, many spectra were obtained at negative epochs. 
For such observations we use a logarithmic-like notation for 
phases, such that 0a = 1.23 means phase -1-0.23 in cycle —1. 

3.2 Ha Properties 

Figure [l] shows the Ha equivalent-width measurements, 
which vary between +1.4 and — 4.3A, folded on the adopted 
ephemeris; it illustrates a number of noteworthy points: 

(i) The Ha light-curve is remarkably symmetrical about 
phase zero. 

(ii) There is a well-defined interval of apparent 'quies- 
cence', lasting ~ 0.3(= l-20o) of the period. 

(iii) The behaviour is repeatable from the earliest quanti- 
tative data (1982)Q supporting a truly periodic underlying 
538-d 'clock'. 

(iv) None the less, the Ha profiles are evidently not 
strictly repeatable; the standard deviation of W\ measure- 
ments about the functional fit is only 0.16A, but scrutiny 



^ This difi^ers from the [Walborn et al.| | |2004| ephemeris (which 
was tied to the less well determined minimum in the Hipparcos 
photometry) by half a cycle. 

^ We measured the equivalent width from a digitized version of 
the plot of Ho given by |Peppel| | |1984[ l, who observed on JD 
2,445,211, 24 years before our last observation. We included this 
point in the fit of eqn. ^ but its exclusion makes no important 
changes to the ephemeris. 
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Figure 1. Ha measurements folded on the ephemeris of eqtn.[2] shown over two cycles. Upper panel: equivalent width (the solid line is 
the ad hoc functional fit, eqtn.[l] described in Section |3.1[ l. Middle panel: FWHM (upper groups of points) and central velocity (lower) 
of excess emission. Bottom panel: the Hipparcos photometry (with a scaled, vertically shifted version of the Ha functional fit to guide 
the eye). 



of the spectra shows that real, small-amplitude variabihty 
contributes to this dispersion at all sufficiently well-sampled 
phases, including quiescence, on timescales longer than a few 
days. We suspect that any 0-type star observed as exten- 
sively and intensively as HD 191612 may show Ha 'jitter' 
at a similar level (cf. 



e.g., 



20041 



Kaper et al. 1997 Morel et al 



The Ha profile during quiescence is asymmetric and filled 
in by stellar- wind emission (Fig. |2|. It is intermediate 
in appearance between the profiles of HD 36861 [A Ori; 
08 ni((f)), pure absorption] and HD 175754 [OS Il((f)), 
P Cygni] , although the closest matches among the spectra at 
our disposal are with HD 193514 [07 Ib(f)] and HD 209975 
[09.5 lb; fig. [2] Thus while we have no exact match to the 
quiescent spectrum, the Ha profile in that state seems to 
be largely unremarkable, with no evidence of substantial ex- 
cess emission relative to normal stars of broadly comparable 
spectral type. 

At other phases, the line-profile morphology is P-Cygni- 
hke, but the increase in emission is not necessarily associated 



® We also examined the spectrum of HD 225160; surprisingly, 
this 08 Ib(f) star shows much stronger emission than the 07 and 
09. 5 lb stars. 



with an increase in the global mass-loss rate (cf. Sec. 4.2 1. 
Phenomenologically, the changes in the appearance of the 
profile can be entirely accounted for by variable amounts 
of roughly gaussian emission superimposed on a constant, 
underlying quiescent-state spectrum (Fig.|2|. We have char- 
acterized the mean velocity displacement and width of this 
excess emission by gaussian fits; results are incorporated 
in Fig. [1] Significant phase dependence is evident for nei- 
ther mean velocity nor fwhm, which average +7 km s~^ 
and 271 km s~^, respectively (standard deviations 10 and 
16 km s~^, commensurate with likely observational uncer- 
tainties); the width of the excess emission is much less than 



both the stellar- wind terminal velocity (Sec. 6.41, and the 
'windy' Ha emission normally seen in luminous O stars. 



Finally, although the Ha and Hipparcos light-curves are 
in phase (to within the errors) , we note that the photometric 
variations are too large to result from excess line emission 
alone; most of the ~3-4% fiux change must arise from true 
continuum-level variations. 
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Figure 2. Balmer-linc profiles for HD 191612 at quiescence 
(0 = 3.61, in black) and near emission maximum {(j> = 0.89, red); 
tickmarks on the y axis are separated by half the continuum level. 
The Ho spectrum of HD 209975 is included to demonstrate that a 
near match to the minimum-state spectrum can be found among 
normal stars (see Section |3.2| for details) . Difference spectra for 
Ha and H(^ are shown at the bottom of the plot. Velocities are 
heliocentric. 



4 THE 538-D PERIOD: OTHER VARIABLE 
LINES 

4.1 Balmer and He I lines 

Other Balmer lines (to at least Hioj^and the He I lines all 
show large variations, in both strength and velocity, which 
are reproducible on the 538-d period. He I 447lA (Fig. [s]) 
exemplifies the typical behaviour of the He I lines; most 
remain in absorption at all phases, although He I AA5876, 
7065, 7281, as well as Hq and H/3, show strong P-Cygni-like 
emission at maximum. While most lines are in absorption 
during quiescent phases (even Ha, as illustrated with other 
Balmer lines in Fig. [2|, He I 6678A is exceptional in retain- 
ing an emission core (Fig. |4|. 



^ Paschen lines from 3-8 to 3-20 are recorded in emission in the 
CFHT ESPaDOnS spectra; these are also strongly variable, in a 
manner consistent with the 538-d period. 



The basic phenomenology of the variability is, evidently, 
infilling of a near-normal underlying absorption profile by 
slightly (~20-30 km s~^) redshifted, almost symmetrical 
emission (cf. Fig. [2|. This interpretation is consistent both 
with direct inspection of line profiles, and with the phase 
dependence of variability; all lines share a common velocity 
around (fja — 0.5, becoming increasingly blueshifted with de- 
creasing line strength (i.e., increasing emission infill). More- 
over, the amplitude of radial-velocity variation correlates 
with the amplitude of line-strength variation. 

This phenomenology appears to be applicable to almost 
all variable features, embracing not only the hydrogen and 
He I lines (including A6678), but also the signature Of?p 
C III 4650A transitions (Fig. |4] the adjacent N in lines are 
much less, if at all, variable, and so are evidently domi- 
nated by 'normal', photospheric, Of emission). He ll A4686 
shows an exceptional pattern of variability; narrow emis- 
sion is present even at quiescent phases, and while the over- 
all emission-line strength increases near phase zero, only 
small changes occur at the position of this narrow emission 
(Fig. [4|, possibly because the emission is already optically 
thick here. 



4.2 UV P-Cygni profiles 



[Walborn et al.| ( j2003[ ) presented a high-resolution Inter- 
national Ultraviolet Explorer (lUE) spectrum obtained on 
1992 Dec 19 (JD 2,448,975.9, 0^ = 9.74, R ~ 10'', AA ~ 
1200-1900 A). There are no other high-dispersion UV spec- 
tra available, but the lUE archive contains a low-resolution, 
short-wavelength spectrum (1984 July 25, JD 2,445,907.0, 
(f)a = 14.04, R ~ 300). After smoothing and binning the 
high-resolution spectrum to render it comparable to the low- 
resolution one (Fig.[5|, there is no evidence for large changes 
in the resonance-line profiles of N v A1240, Si iv A1400 and 
C IV A1550 (although we cannot rule out variations at the 
moderate level observed in the oblique rotator Ori C; 
[Walborn fc Nichols||199^ . This encourages the view that 
the Ha variability results from the changing visibility of a 
constrained plasma, rather than a large-scale global change 
in outfiow characteristics - e.g., if the variability were a con- 
sequence of changes in mass-loss rate, this would be expected 
to have a clear signature in the UV P-Cygni profiles. 



5 THE CONSTANT LINES 

In contrast to the hydrogen and He I lines, the absorption 
lines of metals and of He ll, together with many selective 
emission lines, show only small changes, at most, in line 
strength. For simplicity, we label such lines as 'constant'; if 
there is any line-strength variability, it is at a very low level. 



5.1 Absorption lines 

The C IV 5801, 5812 A doublet exemplifies the behaviour of 
the constant absorption lines. This doublet is particularly 
well-suited to measurement both on astrophysical grounds 
(the lines are very symmetrical, and are formed deep in 
the atmosphere, and hence are less likely to be contami- 
nated by 'windy' emission than many other transitions; cf., 
e.g., Fullerton et al.|["l996 l, and observationally (the nearby 
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Figure 3. Equivalent-width and velocity variations for He I 447lA and He II 454lA. Points with large formal measurement errors are 
omitted for clarity (so that fewer velocities than line strengths are plotted). As in Fig.^ the equivalent- width axis is inverted (such that 
lower points indicate greater absorption/less emission). Points in the left-hand panels are colour coded according to source: E(lodie; |Na2e| 
|et al.|[2007^ , H(igh-resolution), or I(ntermcdiate-dispcrsion). Points in the right-hand panels are colour-coded according to cycle in the 
ephemeris of eqtn. |2] 
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Figure 4. Emission-line variability in HD 191612; spectra are 
labelled by phase in the H« ephcmeris of eqtn.[2] Data obtained 
near Ha maximum are shown in red, those near minimum in 
black; ratio spectra are shown offset by +0.2. (All spectra were 
obtained near orbital zero; Sec. |5.3| ) 

Upper panel: the N III 4640A, C III 4650A, He II 4686A emission- 
line complex; unlabelled tick marks indicate wavelengths of O II 
lines. Note also the infilling of the red wing of He I 4713A. 
Lower panel: He l/He II 6678/6681A. Although A6678 is unusual 
among the He I lines in showing emission throughout 'quies- 
cence', the qualitative nature of the variability (i.e., growth of 
a slightly redshifted, fairly narrow emission superimposed on the 
quiescent-state spectrum) is notable only for the unusually large 
relative amplitude. The AA6717/6722/6729 lines (the latter prob- 
ably a blend in our data, A6728.8-(-6731.1:) are almost constant 
in strength; the positions appear not to be consistent with a stan- 
dard C III identification. 



difltuse interstellar bands at 5778/5780/5797A are of sim- 
ilar strength to the C iv lines, and provide a very use- 
ful zero-point calibration for the wavelength scale, allowing 
rather precise difTerential velocities to be obtained even from 
intermediate-dispersion data of unexceptional quality). 

The C IV measurements are presented in Fig. |6] and 
illustrate the typical behaviour of essentially constant 
line strength coupled with small-amplitude radial-velocity 
variations. The steady increase in velocity between JDs 



2,453,500-700 discussed by Naze et al. (20071 clearly does 



not repeat on the 538-d period. Most metal absorption lines 
strong enough to be measured consistently in most spec- 
tra (e.g., N III AA4511-4534), as well as the He ll absorption 




laOO KiOO 1400 1500 1600 1700 IHOO 

Wavelength (A) 

Figure 5. lUE spectra of HD 191612. Lower spectrum, low- 
resolution, <j>a = 14.04; upper spectrum (offset by 0.5 contin- 
uum units), high-resolution spectrum, (f>a = 9.74, smoothed and 
binned to match the low-resolution data. 



lines (e.g., AA4200, 4541, 5411) follow essentially the same 
behaviour as the archetypal C iv lines, as illustrated in Fig. [3] 
by data for He ll A4541. 



5.2 Selective emission lines 

The high-quality UES and CFHT spectra, in particular, re- 
veal a rich spectrum of selective emission lines (cf. |Walborn| 
20011. In general, the weakness of the lines precludes de- 



tailed scrutiny of their phase dependence in the remaining 
data, but comparison of quiescent and emission-line phases 
in these echellograms establishes that most of the emission 
feature^ exhibit essentially the same behaviour as the C iv 
lines (i.e., near-constant line strength and small-amplitude 
radial- velocity variations). Only the stronger lines can be 
consistently measured, but results for Si iv 6667A shown in 
in Fig.[6]illustrate the general accord between the behaviours 
of the selective emissions and the 'constant' absorption lines. 

The extensive far-red coverage of the high-quality 
CFHT ESPaDOnS spectra allows many other weak emission 
lines to be recognized; these appear also to show little or no 
variability in line strength (at least, between the two epochs 
sampled, at 538-d phases 0.24 and 0.89; the 'variable' lines 
discussed in Section [4] change substantially between these 
epochs). Most of these features are previously unreported 
in 0-star spectra, and currently lack persuasive identifica- 
tions. Measured wavelengths in A (and possible identifica- 
tions, with multiplet numbers from Moore|1945 1 are, for the 
stronger features: 5739.8; 6394.8, 6467.1, 6478.7 (N ill 14? 
good wavelength matches but not all multiplet members 
present); 6482.3; 7002.8; 7037.2 (C ill 6.01); 7306.9; 7455.2; 
7515.7; 8019.2 (N iii 26); 8103.0 (Si iii 37); 8196.6 (C iii 43); 
8251.0; 8265.7; 8268.9; 8286.8; 9705.5, 9715.4 (both strong, 
broad; C ill 2.01). 



Including C II 6578, 6583A; C III 5696A; N II 5001/5005, 5667- 
5680-5686, 6482, 66IOA; N in 5321, 5327A; Si lll(?) 4905. 6A; 
Si IV 6667, 670lA; S IV 4486, 4504A; and the unidentified emis- 
sion lines first noted by[Underhill {199^ at 6717.5, 6722.2, 6729.1 
(probably a blend in our data, 6728.8-f6731.1:), 6744.4A. 
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Figure 6. Equivalent-width and velocity variations for the C IV 580lA absorption and Si IV 6667A emission lines; other details are as 
for Fig|3] 
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Figure 7. Radial-velocity measurements for C IV 580lA (large 
circles), Si IV 6667A (diamonds), and O II lines (small yellow cir- 
cles), plotted with the orbital solution of Table [5] The C IV and 
Si IV points are colour coded according to cycle count in the 
ephemeris of Table [2] and for display purposes the O II veloci- 
ties have been adjusted by —14.3 km s"'^ to bring them to the 
same 7 velocity as the primary (see Sec. |5.3[l. 
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Figure 8. Signatures of the secondary's spectrum. Data have 
been lightly smoothed, rebinned, and renormalized to facili- 
tate comparison between observations, which are labelled by 
orbital/Ho phases (Sec. |5.3[)3.1[ ). Spectra have been shifted to 
bring the primary to zero velocity; long tickmarks indicate rest 
wavelengths of the identified lines. In this reference frame the 
He II, N III and S IV lines show no significant movement (i.e., arise 
principally in the primary's spectrum). Short tickmarks show the 
expected positions of the O II and Si III lines in the secondary 
spectrum, according to the orbital solution given in Table [2] 

5.3 Radial- velocity variations & spectroscopic 
orbit 

The 'constant' lines show significant, systematic radial- 
velocity variationa"] which do not repeat with the 538-d pe- 
riod (Figs. [3] and [6|. These velocity variations cannot reflect 
photospheric motion about a static centre of mass (the im- 

* The narrow core of He II A4686 emission, present throughout 
quiescence, tracks the motion of the 'constant' absorption lines 
and selective emission lines. 



Table 2. Orbital solution. The main orbital parameters are con- 
strained by measurements of C IV 580lA, Si IV 6667A in the pri- 
mary spectrum; K2 is established from O II lines in the secondary 
spectrum. 
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plied change in radius is as great as ~250Rq, or ~17i?*). In 
principle, they could be attributed to changes in the depth 
of line formation in the accelerating outflow, through den- 
sity changes in or near the trans-sonic region resulting from 
stochastic changes in the stellar-wind mass loss rate. How- 
ever, the lack of changes in Ha in the mean quiescent spectra 
around phases 1.5 and 0.5 (corresponding to negative- and 
positive- velocity C iv states) does not encourage confldence 
in this interpretation, and we have already argued that the 
C IV lines are formed relatively deep in the subsonic atmo- 
sphere. 

We therefore consider the possibility of orbital motion 
on periods other than that of the major spectroscopic vari- 
ability. We concentrate on measurements of the C iv 5800A 
doublet and the ~6700A emission-line complex (Si iv + 
unidentifled) , which yield mutually consistent results. Both 
features have velocity zero-points tied in to nearby, mod- 
erately strong diffuse interstellar bands (we adopted DIB 
wavelengths as measured in the CFHT spectra), and there- 
fore record the velocity variations more precisely and accu- 
rately than other lines that lack this advantage. Results are 
given in Table \A2\ 

Fortuitously, the dataset includes good-quality obser- 
vations showing that relatively large positive velocities oc- 
curred in 538-d cycles —3 and —6, as well as at (j>a ~ 0.5, 
hinting at Porb — 3Pa. Trial orbital solutions with Porb — 
1500-1600d yielded residuals that are satisfactorily small, 
but none the less slightly larger than the formal errors re- 
turned by the gaussian fits to lines used to measure ve- 
locities. Since external errors evidently dominate (C iv) or 
match (Si iv) the formal uncertainties, we weighted all C iv 
measurements and all Si iv measurements equally, but with 
the Si IV measurements assigned ~l/3 weight to reflect their 
greater residuals (and adjusted by —2.1 km to bring 
them to the same 7 velocity as the C iv measurements). 

The resulting orbital solution is summarized in Table [2j 
and illustrated in Fig.[7] While this solution has not yet been 
subject to the acid test of predictive power, its success in 
reproducing observations at three separate periastron pas- 
sages, with small residuals that are consistent with realistic 
observational uncertainties, encourages the view that the so- 
lution does characterize a true binary orbit. 
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5.3.1 Secondary orbit 

Careful examination of weak lines of low ionization stages 
in the best data shows velocity displacements in antiphase 
with the stronger lines, bolstering the interpretation of the 
1540-d velocity variations as being orbital in origin. The 
spectroscopic fingerprint of the secondary is weak, but its 
velocity can be reasonably well quantified by simultaneous 
gaussian fits to a selection of a half-dozen unblended, rela- 
tively strong O 11 absorption lines (AA4300-4700A, central 
depths 2-3% below continuum), constrained to share the 
same velocity shift; results are included in Table |A2] 

Although we don't have the benefit of well-defined in- 
terstellar features to provide a velocity reference for most 
observations encompassing the O ll-line spectral region, we 
can exploit the fact that the centrally placed He ll 454lA 
line follows the motion of the C iv and Si iv features used to 



determine the primary's orbit (Sec. 5.1 1. By measuring O ll 



velocity differences with respect to A4541, and correcting for 
computed motion of the primary, we can in effect remove 
spectrum-to-spectrum velocity zero-point offsets. (Using di- 
rectly observed velocities gives essentially identical results, 
but with somewhat larger scatter.) 

An orbital solution of the resulting O ll velocities gives 7 
and K2 (with all other parameters fixed at the primary-orbit 
valueQ; the slope of He 11 vs. He 11— O 11 velocities gives an 
entirely consistent mass ratio. Results are incorporated into 
Table il 



5.3.2 Relationship between the 538- and 1540-d periods? 

The orbital period is ~ 3 x 538d, which has suggested to col- 
leagues the possibility of some sort of resonance. However, a 
solution with Porb = 3Pa gives a significantly poorer fit than 
one with Porh allowed as a free parameter Fj Thus although 
it seems to be quite securely established that HD 191612 is a 
long-period binary, it appears that, with Porb = 2.87Pa, the 
binary orbit has no important role in the large-amplitude 
spectroscopic variability. 

There is, however, an interplay between the periods in 
a limited observational sense, explaining otherwise complex 
behaviour seen in some lines. Fig. [S] shows spectra obtained 
near the extremes of both the orbital motion and the spec- 
troscopic variability. In addition to illustrating the small 
orbital velocity amplitudes (less than the line widths) the 
quiescent-state spectra clearly indicate that absorption in 
the Si III AA4552, 4568, 4575 triplet tracks the secondary 
spectrum, on the orbital period. The primary spectrum has 
emission in these lines which varies on the 538-d period, and 



^ For completeness we report that 7 velocity for the O n lines is, 
formally, -f9.1itl.O km . Although small differences in 7 veloc- 
ities for different lines would not be unusual in an O-type star, it 
should be noted that our 7 velocities for the C IV doublet and the 
Si IV/A6700 complex depend on measured, (not laboratory) wave- 
lengths for both the reference DIBs and the unidentified A6700 
emissions, while the weakness of the O II lines means that small 
zero-point shifts arising from, e.g., unrecognised blends would not 
be surprising. 

An F test yields a probability of <^1% that for the fixed- 
period fit is no poorer than that for the solution with Porb free; 
i.e., formally rules out Porb = 3Pa with > 99% confidence. 




Figure 9. Constraints on the system mass resulting from the or- 
bital solution of Table [2] Solid lines are loci of constant secondary 
mass, labelled in solar masses, implied by the mass function (as- 
suming M2 < J\/i); the dotted line shows the estimated mass 
ratio. 



which shows a strong decline in strength going from A4552 
to A4575 (as is observed in some LBV/WNll spectra; e.g.. 
He 3-519, [Walborn fc Fitzpatrick|2000"| . 

The interplay of these two cycles means that the lines 
can appear purely in absorption (throughout quiescent 538-d 
phases), but in other states the behaviour is more complex. 
When the primary and secondary spectra are at or near 
maximum velocity separation, a P-Cyg-hke appearance re- 
sults in AA4552, 4568, with absorption in A4575 (as in the 
first [top] spectrum in Fig. [sj; emission phases at smaller 
velocity separations can result in apparent disappearance, 
or nulling, of A4552 (third spectrum in Fig. |8|. 

A similar effect is detectable in the N 11 spectrum, in 
particular the 5001/5005A lines, where primary emission 
and secondary absorption clearly move in antiphase; and, 
possibly, in the AA4414.9, 4417.0 O 11 lines. 



5.4 System characteristics 

Mass constraints implied by the orbital solution are illus- 
trated in Fig. |9] We have no constraint on the orbital incli- 
nation, but the masses are entirely consistent with a system 
comprising a ~30-Mq, late-O giant accompanied by a ~15- 
M0, early-B main-sequence star, viewed at an intermediate 
angle (sini ~ 0.5). The projected centres-of-mass separation 
at periastron is (1 — e)asini = (555 ± 68)R0 (~ 40P,); al- 
though this certainly allows for the possibility of wind-wind 
interactions, the Ha emission region, for example, is surely 
much closer to the primary, so that it is unlikely that the 
secondary plays any important role in the 538-d changes. 

The strong variability of the He I classification lines, to- 
gether with the small orbital velocity amplitudes (which en- 
sure that the components' helium lines are never resolved in 
the spectrum) means we cannot estimate the spectral types 
separately from the spectrum; the maximum velocity sepa- 
ration is less than the line widths in the primary. Thus even 
if the components were as just postulated, with an expected 
P-band brightness ratio of ~10:1, the combination of the 
likely differences in spectral types and the resolution and 
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signal:noise of most of our spectra means it is not particu- 
larly surprising that we don't see truly double-lined features 



(other than the peculiar Si in lines noted in section 5.3.21. 
Nonetheless the O ll/Si ill lines indicate a secondary spectral 
type not earlier than BO, and not later than B2 (assuming 
a dwarf luminosity class and normal composition). 



6 OTHER SPECTROSCOPIC PROPERTIES 
6.1 Spectral classification 

Our extended dataset, and improved understanding of the 
spectroscopic variability, allow us to revisit the spectral clas- 
sification (that is, the classification of the spectrum - not the 
component stars). 

The (l>a ~ spectrum (specifically, the WHT spectrum 
of JD 2453941.5) yields spectral type 06.5f?pe ('06.5' from 
the He I 4471/Hell 4541 ratio, with qualifiers 'f?p' from 
C III A4650 and 'e' from the H7 emission). Our refined classi- 
fication of the quiescent-state spectra is 08fp (the weakness 
of C III in that state excluding Of?p). Of course, it is now 
clear that these apparent variability in spectral type does 
not reflect changes in any fundamental stellar properties - 
which would result in large-amplitude photometric variabil- 
ity - but instead results from infilling of the He I A4471 
classification line by emission. 

These results are in full accord with previous classifica- 
tions, including in particular | WalboFni s original blue-region 
spectrum, obtained on 1972 Sept. 24 {<j)a = 22.02) and also 
classified 06.5f?pe (Walborn 19731 - consistent with the 



underlying 'clock' keeping good time for at least a decade 
before the earliest digital data. 

Luminosity classification at any state is not possible 
because of the unique, peculiar profile of the fundamen- 
tal He II 4686A criterion. The strength of Si iv 4089A rel- 
ative to nearby He I lines would be consistent with lumi- 
nosity class V, while the UV Si iv resonance doublet is nei- 
ther dwarf- nor supergiant-like, but is possibly consistent 
with class III. The physical properties established in Sec- 
tion |6]3] are also broadly giant-like. However, all the known 
Of?p stars are sufficiently spectroscopically peculiar that 
any luminosity-class argument based on direct comparisons 
with normal stars must be considered hazardous. 



6.2 Spectroscopic fine analysis 

Although the peculiarities of the spectra preclude extremely 
detailed modelling, the basic photospheric parameters can 
none the less be quite well constrained. We searched for 
matches to the 'quiescent' spectrum of 2 001 Aug 5 



3.61), using a grid of fastwind models (Puis et al. 20051 



sampled at steps of IkK and 0.2 in Tcff and log*; (with Y = 
0.1 and vturh ~ 15 km s~^). We find a generally good match 
with observations for a model with Tcff = 35.0kK, logg = 
3.5, in excellent agreement with the analysis reported by 



Walborn et al. ( 2003 conducted by AH but based on simpler 



models and poorer-quality data). 

No model reproduces the strength of He i 447lA (cf. 
the general discussion of this issue by Repolust et al.|20"04 1 , 
but the He l singlet lines are fully consistent with the esti- 
mated Toff; adopting 34kK would still require logg = 3.5, 




3.5 4,0 4,5 

Log Wavelength (A) 

Figure 10. Spectral-energy distribution for HD 191612 (red). 
The weighted sum of an Ostar2002 model with T^ff = 35kK, 
logg = 3.5 and an Atlas9 model with T^ff = 20kK, logg = 4.5, 
with a V-band flux ratio of 9:1 and E{B — V) = 0.56, is shown 
for comparison (black; see Section |6.3| for details). 



but leads to the model He 11 lines becoming too weak. The 
best-fit model parameters can therefore be considered as 
fairly well determined, to ~ ±lkK, 0.1 in T^ff, logy (al- 
though, of course, no model matches the peculiar He 11 4686, 
He I 6678 A profiles). 

This takes no account of the effect of the secondary 
spectrum, which is expected to have somewhat stronger He I 
lines, while the He 11 line strengths may be underestimated 
by ~10% through dilution. This would result in the primary 
being, if anything, modestly warmer (by perhaps ~lkK) 
than the foregoing analysis suggests. 



6.3 Reddening, radius 

To determine the reddening and angular diameter we com- 
pared a Teft = 35.0kK, loggr = 3.5 Ostar2002 model (|Lanz 



fc Hubenyl |2003[ ) with archival low-resolution lUE spec- 
trophotometry and optical & 2mass photometry. We estl 
mate (E^ff /R0)/(D /kpc) = 6.7 ± 0.05 and E{B - V) = 
0.56 ± 0.03 [using a lCardeUTet al. 1989^ reddening law with 
Rv = A{V)/E(B - V) = 3.1]; here ReS is an 'effective' 
radius characterizing the emitting surfaces. The match be- 
tween the model and observations, though not perfect, is 
reasonably good, and there is no evidence of any IR excess 
out to 2/^m. The fit is slightly improved by allowing for the 
contribution of a cooler secondary, such as a ~20kK sec- 
ondary contributing ~10% of the light at V, and this ad hoc 
model is illustrated in Fig. |10| 

The reddening and general properties are consistent 
with membership of the Cyg OB3 association, as noted by 
Humphreys (19781; we adopt her association distance of 



2.29 kpc (see discussion in | Walborn|2002] for the embedded 
cluster NGC 6871), whence the primary's radius is 



R, ~ 14.5i 



-t-0.2[ 



Ry E(B-V) ] 



JR, 



where /i is the fractional V^-band luminosity of the primary 
and the Teff term accounts for the approximate T^g scaling 
of model-atmosphere V-band surface fiuxes in this temper- 
ature range. The primary then has log(L/LQ) ~ 5.4 and 
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M{V) ~ —5.6 - parameters broadly consistent with a late- 
O giant. 



the effectiveness of the gaussian turbulence in removing in- 
formation with high spatial frequencies from the profiles. 



6.4 Ha mass-loss rate 

In our FASTWIND models the mass-loss rate is characterized 
by a parameter 



Moyr- 



\km [ 



Rr. 



(Puis et al. 



20051, where /d = (p^) / (p)^ is the clumping 
factor; from modelling of Ha, we find — logQ = 12.6-12.7 
for the minimum-state spectrum. We estimate Vao — 
2700 km s"^ (with ~10-15% uncertainty) from the C iv 
resonance doublet in the lUE spectra, whence the primary's 
quiescent stellar-wind mass loss rate is log (^M^ffci^ ~ 

—5.8 dex Meyr"^. As for other aspects of the fastwind 
analysis, this is in satisfactory agreeme nt wi t h the corre- 
sponding result reported by jWalborn et al.| ( 2003 M = 
-5.6 dex Moyr"^). 



6.5 Rotation 

The equatorial velocity expected for a 538-d rotation pe- 
riod and Rt ~ 14.5R0 is ~1.4 km s~^, and a measure- 
ment of Vb sin i would provide a strong test of rotational 
modulation. (For reference, the equatorial rotation velocity 
for synchronous rotation is ~0.5 km s~^, or, for pseudosyn- 
chronous [periastron-synchronized] rotation, ~0.8 km s"'^.) 
Most lines in the spectrum show asymmetries or variability 
to some extent, but, as already discussed, the C iv 5800A 
doublet shows exceptionally symmetrical profiles with little 
or no evidence of contamination by emission; these high- 
excitation lines are expected to form rather deep in the at- 
mosphere, and therefore to be more representative of the 
subsonic photosphere than many other features. 



Donati et al. ( 2006 1; see also Howarth 2003 1 showed 



that the C iv profiles deviate strongly from those expected 
from rotational broadening, but are well matched by sim- 
ple gaussian 'turbulence', with zero rotation. We have at- 
tempted to set upper limits to the rotation rate by look- 
ing for zero-amplitude nodes in the fourier transform of the 



profiles ( |Gray||1992| | in the combined 2005 CFHT spectra 
(s:n > 10'^). We compared the results with transforms of syn- 
thetic spectra that include rotational broadening, isotropic 
gaussian macroturbulence, and gaussian noise, using an in- 



trinsic profile from an Ostar2002 model (Lanz & Hubeny 



20031 



Illustrative results are shown in Fig. |11| Even for these 
very high-quality results there is no sign of rotational ze- 
roes above the noise (which dominates at frequencies > 
10~^ s km~^). Empirically, we find that the observations 
are well matched by any model having 



\J(^Y + {vc sini/2)2 



45 km s 



(where cry characterizes the velocity dispersion of the macro- 
turbulence and Vc sin i characterizes the projected equatorial 
rotation velocity), for WeSini<60 km s~^; larger values of 
Vc sin i are ruled out by the data. This is a disappointingly 
weak upper limit given the data quality, but results from 



7 DISCUSSION 

The symmetry and reproducibility of, in particular, the Ha 
538-d light-curve strongly suggest an origin in an essentially 
geometrical process - i.e., the changing aspect of a distinct 
emission region. This is also consistent with the absence of 
evidence of associated global changes in stellar-wind proper- 



ties (Sec. 4.2 1. Since we have found that the variability is not 
orbital, our results lend strong support to the proposal by 



Donati et al. ( 2006 1) of rotationally modulated line emission 
from a magnetically constrained plasma. 

One might hope that the fairly rich emission spectrum 
would offer insight as to conditions in the line-forming re- 
gion, but there is a dearth of traditional nebular diagnostics 
(suggesting a high-density regime). We can only infer some 
general characteristics of the Ha-emitting region. First, op- 
tically thin Ha emission in an any time-independent ax- 
isymmetric structure cannot account for the observed 538-d 
behaviour; at least half the emission from any such struc- 
ture would visible at all times, but emission equivalent to 
the adopted 'excess' emission cannot plausibly be present 
during 'quiescent' phases. At the same time, very complex 
geometries (such as observed in r Sco; Donati et al.|2006 3) 
are unlikely, as they would not be consistent with the large 
amplitude of emission variability. 

Secondly, for the distance and reddening adopted in Sec- 
tion [6]3j and assuming isotropic, case-B recombination emis- 
sion, the maximum excess Ha emission over the quiescent 
state corresponds to 



nidV : 



7 X lO^^cm" 



This is a lower limit to the true emission measure if the 
emission is not optically thin. 

If we adopt the hypothesis that the geometry of the 
line-emitting region is determined by the magnetic field, 
then the simplest acceptable magnetic-field geometry is a 
misaligned, centred dipole - known to provide a reasonable 
approximation to Ori C, the prototype magnetic O star 



(!Stahl et al.||1996' Donati et al.||2002| iWade et al.||2006l 



|Babel fc Montmerle, (1997 1 argue that, for such a geometry, 
stellar-wind material from the magnetic polar and temper- 
ate regions will be defiected along the field lines towards the 
magnetic equator, where the colliding fiows shock to create a 
high-temperature. X-ray emitting plasma. The plasma cools 
to form a dense disc which gives rise to optical emission. 
MHD calculations support this general outhne, albeit with 
a number of refinements (e.g., |ud-Doula & Owocki||2002 



Gagne et al.|2005[|ud-Doula et al.|2006| ). 

Since this model is at least consistent with the direct 
magnetic-field measurements of HD 191612 ( [Donati et al. 
[2006 a) , for heuristic purposes we explore schematic 'toy' 
models of Ha emission from a centred, tilted, geometrically 
thin disk, in which the relative emission is simply related to 
the projected area (taking into account limb darkening and 
occultation by the star). Although this is a physics-free, ge- 
ometrical model, we can speculate that it may characterize 
a more realistic scenario; e.g., for a disk-like emission region. 
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Figure 11. Illustrative line-broadening models for the C IV 580lA line. 

Upper panels: results in the fourier domain. Solid black lines show the transform of the observed profile (i.e., the transform of the black 
spectrum in the lower panels); red lines show the transforms of synthetic profiles computed with isotropic gaussian macroturbulence 
(with dispersion a, in km s~^), rotational broadening (with iicsini = v), and gaussian noise; and green lines show the transforms for 
rotational broadening alone. Any model with tie sin j< 60 km s~^ (and appropriate macroturbulence) is acceptable. 

Lower panels: results in the wavelength domain. The observed spectrum is shown in each panel in green (the directly observed profile) 
and black (the version of the profile used in the analysis, after re-rectification to take out the diffuse interstellar band). The synthetic 
profiles, shown in red, are the inverse transforms of the models shown in red in the upper panels. At this scale, the synthetic profiles 
(which have been modestly scaled in intensity to match the observed line strength) are almost indistinguishable from each other, and 
from the observed spectrum, in the wavelength domain. 
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Figure 12. 'Toy' models compared to phased Ha equivalent 
widths (circles). Lower model (left-hand axis), a single surface 
'spot'; upper model (right-hand axis), an optically thick disk. See 
section [t] for details. 



the 'hmb darkening' may actually correspond to increasing 
Ha optical depth as the line of sight approaches the plane 
of the disk. 

We find that this simple model can provide a reasonable 
match to the Ha variability, provided that (i) the sum of the 
inclination of the rotational axis to the line of sight and the 



angle between that axis and the magnetic axis is close to 
90°, and (ii) moderately strong limb-darkening is present 
(so that the emission is low at all rotational phases when 
the line of sight is close to the disk plane). Moreover, the 
width:thickness ratio must be reasonably large, to account 
for the large amplitude of emission-line variability. 

A specific model, selected through a genetic-algorithm 
minimization, is shown in Fig. |12| its parameters are in- 
ner radius i?in = 2.00i?, , outer radius i?out = 2.13i?*, axial 
inclination i — 68°, magnetic-axis offset a = 27°, and lin- 
ear limb-darkening coefficient u = 0.7. These parameters 
should be regarded only as illustrative, as many other com- 
binations give closely similar fits (e.g., 1.55iZ,, 1.76R*, 41°, 
56°, 0.6 gives results almost indistinguishable from those in 



Fig. 121. The radii, in particular, are only very weakly con- 



strained from the Ha light-curve, because the total emission 
is fixed by an arbitrary scaling factor (parametrizing the sur- 
face emissivity). 

[In principle, the contiimum photometry could help fix 
the radii, but in practice it only weakly constrains the model, 
because the photometric amplitude is so small, the noise is 
relatively high, and the number of free parameters is large. 
The only firm conclusion is that a very large disk, very close 
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to the star, and optically thick in the continuum, is not 
allowed. The broad-band photometry suggests that any disk 
cannot contribute more than ~3-4% of the visible-region 
broad-band (Hp) flux - that is, the implied equivalent width 
of the Ha emission referred to disk continuum is as much as 
~150A.] 

The toy model succeeds in reproducing the Ha light- 
curve, but it would clearly be rash to put too much weight 
on this success. Not only does the physical model which 
inspired it encounter difficulties in accounting for the details 
of the observed X-ray emission ( Naze et al.||2007 l, but also 
the disk model is certainly not unique. To illustrate this, 
we have also considered a minimalist model of Ha emission 
from a single surface spotj^ described by 

fa = fo + ^(cos i cos /3 + sin i sin P cos — u + u cos /i), 

where fo,A are normalizing constants, i is the inclination of 
the rotation axis to the line of sight, (3 is the colatitude of the 
spot, 4> is the rotational phase, m is a linear limb-darkening 
coefficient, and fi is the angle between the line of sight and 
the surface normal at the spot. We find that the form of the 
Ha light-curve can again be well matched by this model, 
with u ~ 0.6, i + /3~ 105°, i ~ f3 (Fig.[l2|. 



A continuing programme of circular spectropolarimetry is 
working towards this aim, and a magnetic-field measure- 
ment from the second-epoch ESPaDOnS observation used 
here already gives a clearly detected field with a longitudi- 
nal component larger than the discovery measurement, as 



expected for the geometry adopted by Donati et al. ( 2006 1) 
and in our toy disk model. 



ACKNOWLEDGEMENTS 

We thank all the observers who contributed spectra to this 
campaign; the service observers at the Isaac Newton Group; 
Sergio Ilovaisky and Ashok K. Pati for orchestrating spec- 
troscopy at OHP and VBO; Roger Griffin for comments on 
the orbital solutions; and our referee, Otmar Stahl, and his 
colleagues at the Landessternwarte Heidelberg for establish- 
ing the date of Peppel's observation from the original logs. 
We also gratefully acknowledge the following for support: 
the FNRS, PRODEX/Belspo, and OPTICON (YN, GR); 
the Estonian Science Foundation, Grant 6810 (KA, IK); and 
the F.H. Levinson fund of the Peninsula Community Foun- 
dation (DM). 



8 SUMMARY AND CONCLUSIONS 
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APPENDIX A: OBSERVATIONAL DETAILS 

This appendix summarizes 



line) 



The log of observations (Table Al available in full on- 

) 

Radial-velocity measurements used in the orbital solu- 



tions (Table A2 1 
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Table Al. Format of Table Al (available in full on-line), the log of optical spectroscopy used in this paper. Each 'observation' is a 
datasot collected at a single site in a given night, and may consist of several separate spectra. Phases are calculated with the ephemerides 
of TablejSjand eqtn.jS] the notation C.ppp means phase +.ppp in cycle —C. Where two observers are identified, the second name is the 
PI or instigator of the observation; 'srv' indicates a service-mode observation. Where observers provided smoothed or binned spectra, 
the effective resolution is tabulated (AA). 



JD 


Year 


Phase (pa 


Phase (fi 




A range 




Telescope / 


Observer 


AA 








(Ha) 


(orb) 




(nm) 




Instrument 




(A) 


(A) 


2447691.5 


1989.45 


TT.35 


4.09 




651- 


-662 


INT/IDS 


Prinja/Howarth 


0.5 


-1-1.37 


2447724.7 


1989.54 


TT.42 


4.11 


404- 


-499 




INT/IDS 


Herrero 


0.8 




2447726.5 


1989.55 


TT.42 


4.11 




634- 


-677 


INT/IDS 


Herrero 


1.5 


-1-1.41 


2448117.4 


1990.62 


TO.15 


4.37 


394- 


-474 




INT/IDS 


Howarth 


1.5 




2449139.7 


1993.41 


8.05 


3.03 


439- 


-481 




INT/IDS 


Vilchez 


0.5 




2453955.4 


2006.60 


1.00 


0.15 




637- 


-673 


Tartu 


Kolka 


1.5 


-4.00 


2453956.5 


2006.60 


1.01 


0.15 




637- 


-673 


Tartu 


Kolka 


1.5 


-3.96 


2453958.4 


2006.61 


1.01 


0.15 




637- 


-673 


Tartu 


Kolka 


1.5 


-4.14 


2453966.6 


2006.63 


1.03 


0.16 




384-706 




WHT/ISIS 


Leisy/Lennon 


0.4 


-4.07 


2454002.5 


2006.73 


1.09 


0.18 




384-706 




WHT/ISIS 


Leisy/Lennon 


0.4 


-3.32 



Table A2. Velocities used in orbital solutions. Note that each of the three datasets is subject to a separate velocity zero-point error of 
order a few km s~^ (see footnote IM . 



JD 


Vq 


JD 


Vq 


JD 


V© 


JD 


Vq 


JD 


Vq 




km s ^ 




km s ^ 




km s ^ 




km s ^ 




km s ^ 


C IV 5800A 


2449529.4 


-9.3 


2453203.5 


-11.3 


2453275.2 


-11.3 


2453474.6 


-8.8 


2453601.4 


4-4.4 


2452127.5 


4-9.0 


2453225.8 


-6.3 


2453276.2 


-10.6 


2453488.6 


-7.3 


2453625.4 


4-3.7 


2452528.8 


-8.3 


2453246.3 


-8.9 


2453278.3 


-11.7 


2453519.6 


-5.9 


2453652.3 


4-7.0 


2452549.4 


-3.1 


2453251.3 


-9.6 


2453282.4 


-12.0 


2453536.7 


-5.8 


2453653.3 


4-7.6 


2452607.3 


-5.3 


2453264.3 


-7.8 


2453283.3 


-10.1 


2453545.0 


-3.2 


2453668.3 


4-6.3 


2452801.6 


-11.4 


2453265.3 


-9.6 


2453302.3 


-11.3 


2453546.0 


-3.7 


2453681.3 


4-11.7 


2452808.5 


-11.2 


2453266.3 


-12.2 


2453314.6 


-10.3 


2453547.0 


-3.5 


2453836.8 


4-6.3 


2452816.6 


-11.0 


2453267.2 


-9.9 


2453316.3 


-9.7 


2453548.0 


-4.3 


2453892.7 


-0.2 


2452871.2 


-12.1 


2453267.4 


-8.2 


2453320.3 


-12.9 


2453554.6 


-7.4 


2453896.0 


+2.1 


2453132.9 


-6.5 


2453269.2 


-10.9 


2453324.3 


-12.7 


2453568.5 


-0.8 


2453941.5 


-3.4 


2453146.7 


-10.4 


2453270.2 


-14.2 


2453326.3 


-13.2 


2453585.5 


-2.2 


2453966.6 


-1.7 


2453192.5 


-7.6 


2453272.4 


-8.3 


2453346.4 


-9.5 


2453594.5 


-4.5 


2454002.5 


-3.2 


A6700 


emission 


















2450683.5 


-9.4 


2452808.5 


-34.8 


2453266.6 


-36.6 


2453474.6 


-33.3 


2453594.5 


-22.0 


2452127.4 


-16.4 


2452919.3 


-34.6 


2453270.2 


-30.2 


2453480.7 


-29.6 


2453607.4 


-15.6 


2452128.4 


-16.2 


2453132.9 


-42.5 


2453272.4 


-31.4 


2453488.6 


-33.6 


2453625.4 


-24.2 


2452129.4 


-17.6 


2453146.7 


-38.8 


2453278.3 


-37.8 


2453519.6 


-26.6 


2453652.3 


-18.2 


2452129.6 


-15.3 


2453157.7 


-35.3 


2453283.3 


-36.0 


2453536.7 


-16.6 


2453653.3 


-23.1 


2452130.4 


-15.6 


2453209.7 


-37.2 


2453290.6 


-34.7 


2453545.0 


-28.0 


2453668.3 


-16.5 


2452132.7 


-18.3 


2453210.5 


-37.6 


2453291.6 


-36.2 


2453546.0 


-28.0 


2453681.3 


-23.1 


2452483.6 


-29.9 


2453225.8 


-33.2 


2453326.3 


-41.4 


2453547.0 


-27.8 


2453836.8 


-16.1 


2452566.3 


-29.2 


2453246.3 


-40.0 


2453368.3 


-31.4 


2453548.0 


-27.2 


2453892.7 


-27.2 


2452607.3 


-32.0 


2453250.5 


-38.2 


2453369.3 


-37.7 


2453554.6 


-23.3 


2453896.0 


-17.5 


2452801.6 


-32.6 


2453251.3 


-41.5 


2453446.7 


-39.0 


2453585.5 


-27.5 


2453966.6 


-25.8 


O II absorption 


2452101.6 


-28.7 


2452566.3 


-24.6 


2453157.7 


-f21.6 


2453266.6 


4-19.7 


2453681.3 


-35.1 


2452127.4 


-33.0 


2452632.3 


-5.9 


2453209.7 


-fl5.7 


2453554.6 


-6.6 


2453836.8 


-32.7 


2452127.5 


-33.5 


2452808.5 


-1-7.2 


2453210.5 


-f6.9 


2453568.5 


-11.1 


2453892.7 


-21.4 


2452128.4 


-33.4 


2452816.6 


+5.7 


2453225.8 


-fll.3 


2453625.4 


-26.3 


2453896.0 


-23.4 


2452129.4 


-32.7 


2452831.6 


+6.1 


2453246.3 


-1-5.0 


2453652.3 


-19.1 


2453941.5 


-17.2 


2452130.4 


-34.9 


2452837.6 


-1-14.3 


2453250.5 


4-11.6 


2453653.3 


-25.8 


2453966.6 


-20.4 
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Table Al. Optical spectroscopy used in this paper. Each 'observation' is a dataset collected at a single site in a given night, and may 
consist of several separate spectra. Phases are calculated with the ephemerides of Table 2 and eqtn. 2; the notation C.ppp means phase 
+.ppp in cycle — C. Where two observers are identified, the second name is the PI or instigator of the observation; 'srv' indicates a 
service-mode observation. Where observers provided smoothed or binned spectra, the effective resolution is tabulated (AA). 



JD 


Year 


Phase 4>a 
(Ha) 


Phase If 
(orb) 




A range 
(nm) 


Telescope/ 
Instrument 


Observer 


AA 

(A) 


(A) 


2447691.5 


1989.45 


TT.35 


4.09 






651-662 


INT/IDS 


Prinja/Howarth 


0.5 


H-1.37 


2447724.7 


1989.54 


TT.42 


4.11 


404- 


-499 




INT/IDS 


Herrero 


0.8 




2447726.5 


1989.55 


TT.42 


4.11 






634-677 


INT/IDS 


Herrero 


1.5 


+1.41 


2448117.4 


1990.62 


TO.15 


4.37 


394- 


-474 




INT/IDS 


Howarth 


1.5 




2449139.7 


1993.41 


8.05 


3.03 


439- 


-481 




INT/IDS 


Vilchez 


0.5 




2449529.4 


1994.48 


8.77 


3.28 


445- 


-603 




WHT/ISIS 


Crowther 


1.5 




2450683.5 


1997.64 


6.92 


2.03 


356- 


-407 


641-677 


WHT/UES 


Trapero 


0.07, 0.13 


-3.82 


2452101.6 


2001.52 


3.56 


2.95 


369- 


-535 




INT/IDS 


Erwin /Herrero 


0.20 




2452127.4 


2001.59 


3.60 


2.97 


388- 


-482 


636-676 


WHT/lSlS 


Evans /Howarth 


0.5, 0.8 


-1-1.58 


2452127.5 


2001.59 


3.60 


2.97 


361- 


-599 




WHT/UES 


Evans /Howarth 


0.06 




2452128.4 


2001.60 


3.61 


2.97 


388- 


-477 


636-676 


WHT/ISIS 


Evans /Howarth 


0.5, 0.8 


-1-1.59 


2452129.4 


2001.60 


3.61 


2.97 


388 


478 


636-676 


WHT/ISIS 


Evans /Howarth 


0.5, 0.8 


-1-1.51 


2452129.6 


2001.60 


3.61 


2.97 


365- 


-543 


636-676 


WHT/ISIS 


Evans /Howarth 


0.8, 0.8 


-1-1.46 


2452130.4 


2001.60 


3.61 


2.97 


388- 


-478 


636-676 


WHT/ISIS 


Evans /Howarth 


0.5, 0.8 


-1-1.49 


2452132.7 


2001.61 


3.61 


2.97 






635-677 


INT/IDS 


Srv/Herrero 


0.6 


-1-1.35 


2452483.6 


2002.57 


2.27 


T.20 






625-681 


INT/IDS 


Srv/Herrero 


0.6 


-1-0.39 


2452508.2 


2002.64 


2.31 


1.21 




322- 


-703 


INT/IDS 


Arauj o/Gansicke 


3.8 


-1-0.79 


2452528.8 


2002.69 


2.35 


T.23 




320- 


-733 


Bok 


Wagncr/Bond 


1.6, 2.6 


-1-1.19 


2452549.4 


2002.75 


2.39 


1.24 




440- 


-775 


WHT/ISIS 


Rix/Pcttini 


1.8-2.6 


-1-1.34 


2452566.3 


2002.80 


2.42 


T.25 


395- 


-485 


633-708 


WHT/ISIS 


Crowther 


0.8 


-1-1.42 


2452607.3 


2002.91 


2.50 


T.28 




382- 


-797 


WHT/ISIS 


Evans /Crowther 


0.8, 3.0 


+1.33 


2452632.3 


2002.98 


2.54 


T.29 


405- 


-495 




WHT/ISIS 


Harries 


0.7 




2452775.0 


2003.37 


2.81 


1.39 


405- 


-473 




WIYN 


Bond 


0.6 




2452801.6 


2003.44 


2.86 


1.40 




386- 


-684 


WHT/ISIS 


Lennon /Herrero 


0.6-1.2 


-2.70 


2452803.7 


2003.45 


2.86 


1.41 


425- 


-472 




INT/IDS 


Christian 


0.3 




2452808.5 


2003.46 


2.87 


T.41 




380- 


-700 


WHT/ISIS 


Herrero 


0.7 


-2.66 


2452810.8 


2003.47 


2.88 


T.41 


405- 


-473 




WIYN 


Bond 


3.9 




2452816.6 


2003.48 


2.89 


T.41 




330- 


-684 


INT/IDS 


Lennon 


0.6 


-2.67 


2452831.6 


2003.52 


2.91 


T.42 


372- 


-519 


599-697 


INT/IDS 


Prada 


0.7 


-3.00 


2452832.5 


2003.53 


2.92 


T.42 


332- 


-498 


586-705 


INT/IDS 


Howarth 


1.4, 1.1 


-3.69 


2452833.4 


2003.53 


2.92 


T.42 






587-704 


INT/IDS 


Howarth 


1.1 


-3.92 


2452834.7 


2003.53 


2.92 


1.43 


353- 


-521 




INT/IDS 


Howarth 


1.1 




2452835.7 


2003.53 


2.92 


1.43 


353- 


-522 


585-697 


INT/IDS 


Howarth 


1.4, 1.1 


-3.79 


2452836.7 


2003.54 


2.92 


T.43 


352- 


-522 




INT/IDS 


Howarth 


1.4 




2452837.6 


2003.54 


2.93 


T.43 


382- 


-490 


630-683 


INT/IDS 


Howarth 


1.1, 0.5 


-3.74 


2452838.6 


2003.54 


2.93 


T.43 


382- 


-491 


627-683 


INT/IDS 


Howarth 


1.1, 0.5 


-3.15 


2452845.9 


2003.56 


2.94 


T.43 


405- 


-473 




WIYN 


Harmcr/Bond 


0.6 




2452854.8 


2003.59 


2.96 


T.44 


405- 


-473 




WIYN 


Harmer/Bond 


0.6 




2452871.2 


2003.63 


2.99 


T.45 




391- 


-716 


OMM 


Pellerin 


3.9 


-4.30 


2452902.6 


2003.72 


1.05 


T.47 


405- 


-473 




WIYN 


Bond 


0.6 




2452913.7 


2003.75 


T.07 


1.48 


334- 


-648 




MMT 


Bond 


3.6 




2452919.3 


2003.76 


1.08 


T.48 






634-677 


OHP/Aurelie 


Rauw 


0.6 


-3.61 


2452922.4 


2003.77 


T.08 


T.48 


446- 


-490 




OHP/Aurelie 


Rauw 


0.6 




2452980.3 


2003.93 


1.19 


T.52 


373- 


-547 


576-756 


WHT/ISIS 


Lennon /Herrero 


1.2 


-1.40 


2453132.9 


2004.35 


T.47 


T.62 




444- 


-674 


OHP/Elodie 


Rauw 


0.3 


+1.52 


2453132.9 


2004.35 


T.47 


T.62 


405- 


-474 




WIYN 


Harmer/Bond 


0.6 




2453146.7 


2004.39 


T.50 


T.63 


374- 


-550 


568-742 


WHT/ISIS 


Licandro 


1.6 


+1.30 


2453147.4 


2004.39 


T.50 


T.63 




371- 


-693 


Skinakas 


Reig 


2.0 


+1.27 


2453157.7 


2004.42 


T.52 


T.64 


359- 


-494 


604-681 


WHT/ISIS 


Howarth 


0.8 


+1.30 


2453180.4 


2004.48 


T.56 


T.65 






550-686 


Skinakas 


Reig 


2.0 


+1.30 


2453182.4 


2004.48 


T.57 


T.65 


375- 


-510 




Skinakas 


Reig 


2.0 




2453192.5 


2004.51 


1.59 


T.66 




390- 


-682 


OHP/Elodie 


Siviero 


0.1 


+1.29 


2453193.4 


2004.51 


T.59 


T.66 






551-687 


Skinakas 


Reig 


2.0 


+1.45 


2453194.4 


2004.52 


T.59 


T.66 


377- 


-509 




Skinakas 


Reig 


2.0 




2453202.5 


2004.54 


T.60 


T.66 


350- 


-530 


609-818 


Loiano 


Negueruela 


3.0, 3.0 


+1.50 


2453203.5 


2004.54 


T.61 


T.66 




388- 


-999 


Loiano 


Negueruela 


0.7 


+1.47 


2453209.7 


2004.56 


1.62 


1.67 


449- 


-540 


597-688 


WHT/ISIS 


R.C. Smith 


0.7 


+1.36 


2453210.5 


2004.56 


T.62 


T.67 


449- 


-540 


597-689 


WHT/ISIS 


R.C. Smith 


0.7 


+1.34 
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JD 


Year 


Phase (pa 


Phase ifi 




A range 


Telescope/ 


Observer 


AA 


WA(Ha) 






(Ha) 


(orb) 




(nm) 


Instrument 




(A) 


(A) 




9nn4 fin 


1.65 


T.68 




372- 


-684 


SPM 


VjrcOl git; V 


0.4 


+1.31 


2453243.4 


2004 65 


1.68 


1.69 




362- 


-684 




Rpiff 


2.0 


+1.06 




9(104 RR 


1.69 


T.69 


401- 


-590 


653-674 


OHP /Rlnrlip 


01 V j JTvcLU W 


0.3 


+1.02 


2453250.5 


2004 67 


T.69 


1.70 


353- 


-505 


602-688 


WHT/TRTS 


Ostcnscn /Howcirth 


0.6 


+0.89 




9004 R7 
zuu^. u / 


1.70 


T.70 


401- 


-590 


653-674 


OFTP /FlnHip 


01 V / ±V,(1.L1W 


0.3 


+0.89 


94'^'^9'^9 A 


9004 R7 
zuu^. u / 


1.70 


T.70 




365- 


-683 


Iv 1 i i divirl 3 


Rcig 


2.0 


+0.90 


94.'^'^9fil 


9004 70 
ZUU^. 1 u 


1.71 


T.70 


363- 


-570 




Ivl 1 1 a. Jvclo 


Reig 


2.0 




O/II^QOfiO /I 

Z4:ijO.^OZ .'1 


9004 70 
ZUU^. f u 


T 79 


T 7n 




363- 


-684 


Skina.k£LS 


Reig 


9 n 


-UO I'X 


Zi t !j O U . O 


2004.71 


1.72 


T.70 




420- 


-710 


1 .T"! TYl 


r\. llULIivlllll 


0.5 


+0.72 


24'i32fi'i S 


2004.71 


1.72 


T.71 




420- 


-710 


V—'l IIUCCL 


A ntrttrnn 

IX U 1^ J\.XX1 IX 


0.5 


+0.62 


94c;o9cc o 

Z'itJOZUU.O 


2004.71 


1.72 


T.71 




420- 


-710 


1 .Ti im 


A "rri"/"!!^ rii Ti 


0.5 


+0.53 


94'=l'?9fifi R 

ZttJOZrUU .U 


2004.71 


1.72 


T.71 


346- 


-505 


597-689 


WHT /T*^m 

VV n J- / loio 


r^^iTin / rTrf^THra vi" rl 
UCllll / llLIWdl lili 


0.6 


+0.58 




9004 79 


1.72 


T.71 




420- 


-710 


« r 1 m r» 

v.-' 1 iillL^cX 


A Tl + r^V rl 1 11 
Jill LvJivlllll 


0.5 


+0.58 


2/1 ■S'^2fi7 4 


2004 72 


1.73 


T.71 


401- 


-590 


653-674 


OHP/Klodip 


'~lT'\r / R il 1 \WT 
01 V / XVcjjUW 


0.3 


+0.60 




9004 79 


T.73 


T.71 




420- 


-710 


vyl lllltJirl 


A yi'i"*^Lrrl1Tl 
rt.ll LVJlVllill 


0.5 


+0.62 


Z^U OZ 1 u . z 


9004 79 


T 7*^ 


T 71 




420- 


-710 




Antokhin 


u.o 


-UO 40 


94'^'^979 4 


9004 7*^ 


1.73 


T.71 


401- 


-590 


653-674 


OHP /F,lnrlip 


01 V j JvdU-W 


0.3 


+0.54 




2004.74 


T.74 


T.71 




420- 


-710 


1; 1*1X1(1*311 


A Titr»k"niTi 

l\. IX L ^IVXXX IX 


0.5 


+0.30 


2453276 2 


2004.74 


T.74 


T.71 




420- 


-710 


(It"! m pa 


A nt okhi n 


0.5 


+0.18 


94r;Q97Q q 
Z'iOOZ / o.o 


9004 7^ 


T 7^^ 

J.. 1 


T 71 


401- 


-590 


653-674 


OHP /PlrtHiP 

wiii y rLiiociitJ 


S r V / RcLUW 


u.o 


4-0 9^^ 


Z^tJOZOZ.'if 


9004 7R 


1.75 


T.72 




365- 


-685 






2.0 




94'=i'^989 4 


9004 7R 


1.75 


T.72 




370- 


-790 




A Tl von '7'71 


0.04 


+0.22 


ZtiJijZOO .o 


9004 7R 


1.75 


T.72 


401- 


-590 


653-674 


OHP /FloHip 


01 V j JXcl.LlW 


0.3 


+0.18 


24^'^2Sfi 4 


2004.77 


T.76 


T.72 


445- 


-490 






X y.(X Li W 


0.6 




9zL'^'^9^^Q '\ 
z^oozoy . o 


9004 79, 


1.77 


T.72 






636-675 


OHP / Anrplip 

VJlliT / iT. Lll die 


Ra.iiw 


0.6 


—0.27 




9004 78 


T.77 


T.72 






647-714 


KPNO 


Gies 


0.7 


—0.42 


2453291.6 


2004 78 


T.77 


T.72 






647-714 


KPNO 


Gies 


0.7 


—0.31 


24532Q3 3 


2004 79 


1.77 


T.72 


445- 


-490 




OHP/Aiirelie 


X bell u w 


0.6 




94'=l'^^n9 ^ 

Z'ltJOiJUZ . 


9004 81 


1.79 


T.73 


401- 


-590 


653-674 


OHP /F,lnrlip 

KJ 111 / JZiiUUitJ 


01 V / rVcLUW 


0.3 


—0.64 


24533(14 3 


2004 82 


1.79 


T.73 






644-662 


^tlri n a li"ji <3 

Ivl 1 1 Cli2VCU3 




2.0 


—0.80 


94'=l'^m A fi 

Z 'itJOtJ -L ^ • U 


9004 84 


1.81 


T.74 




366- 


-668 


SPM 


Vjrcwl git; V 


0.3 


— 1.44 




9004 8"^ 


1.82 


T.74 


401- 


-590 


653-674 


(-jtrp /plfirlip 
v_/jrLX^ / njiouic 


Oi V / JTvcLUW 


0.3 


— 1.44 


94c;QQ9n s 


9004 Sfi 


1.82 


T.74 


401- 


-590 


653-674 


OHP /RloHip 

V./±l± / l_Jl*_J*_J.lt> 


^T'^r / R a 1 nsr 
1 V / X vt:i Li w 


0.3 


— 1.60 


94c;qq94 q 


9004 R7 


1 .00 


T 7A 




390- 


-682 




S r V / Rciu w 


n 1 
u. 1 


1 77 


z'rtuoozu .0 


9004 9S< 


T.83 


T.74 


401- 


-590 


653-674 


OHP /FloHip 


01 V / J.\i<iu.w 


0.1 


— 1.85 


94'^'^'?4R 4 


9004 Q'^ 


1.87 


T.76 




381- 


-703 


Pii-TV/Tirti /A/Tiicifnc 

JT H..,iVliLtl/ IVl Ll&ll^L'& 


IN eg U-Cl U-Cidi 


0.1 


—3.25 


ZIUOOUO-O 


9004 QQ 


1.91 


T.77 


360- 


-537 


613-704 


WHT /TSTS 
VV n _L / lOJ-O 


T-</^o1j^otc / ^TJ^'^xr^^ vi" ri 
l-VJClCdo / J-lCWdil lill 


1.1, 0.6 


—3.65 




2004 QQ 


1.91 


T.77 


360- 


-537 


613-704 


WHT/TRTS 


Rr»plr^pfc! /H^mxrav+ri 

XX.IJCX'JCXO/ J-X^JWCtX lilX 


1.1, 0.6 


—3.53 


94'=1'^'^71 ^ 

Z^OOO / X .0 


900^ 00 


1.92 


T.77 


353- 


-502 


611-777 


WHT /TSTS 


vjr cLtS 11 b 1 Ulvc 


1.7 3.2 


—3.48 


94'=1'^444 « 


900^ 90 


0.06 


T.82 




364- 


-834 


NOT/ AT,FnS(^ 


XjlUcLllUl U 


2.8 


—3.90 


94'=1'^44R 7 


900"^ 91 

Z\JL>iJ. Z J. 


0.06 


T.82 






639-705 


3 1 ctg U 


Ivlorcl 


1.2 


—3.53 


94'^'^4'^R ^ 


2005 24 


0.08 


T.83 






637-673 


Tartu 


A Ti m 1 It / PC If CI 

XA.1111 LliV / XVWliVtX 


1.5 


—3.70 


94'^'^4R'^ R 


900*^ 9R 
zuuo . zu 


0.09 


T.83 


401- 


-590 


653-674 


QUp /PloHlP 

V_yl.lJr / ll/ltJU-lC 


^"VAT" / R ^1 1 

01 V / IXdUVV 


0.3 


—3.22 


24'^ '^479 8 


900^ 9S 


0.11 


T.84 






550-700 


FMO 


Mr David 


1.8 


—3.29 


2453474 6 


2005 28 


0.11 


1.84 


401- 


-590 


653-674 


OHP /Elodip 


Srv /RsjUw 


0.3 


—3.01 


24534Sn 7 


2005 30 


0.12 


T.84 


350- 


-505 


600-689 


WHT /TSTS 


Licajiidro 


0.6 


— 2.96 


2453486 Q 


2005 32 


0.13 


1.85 






560-680 


FMO 


IVTrDaviH 

xvx ct V xvj 


1.8 


—2.80 


941=1 '^4J<S R 
z^oo^oo.u 


900^ ^9 


0.14 


T.85 


401- 


-590 


653-674 


OHP /FlnHip 


01 V / rXiCLLiW 


0.3 


—2.42 




200^ SQ 


0.19 


1.87 






509-716 


Sllri n a <; 


Rpip- 
rteig 


2.0 


— 1.49 




900^ 41 


0.19 


T.87 


401- 


-590 


653-674 




Oi V / JTvcLUW 


0.3 


—1.08 




900"^ 4"^ 


0.23 


T.88 


401- 


-590 


653-674 


OHP /FlnHip 

V_yl.lJr / l_jlwi_llt^ 


01 V y' JXd.LlW 


0.3 


—0.44 


94'^'^'^4'^ n 


900"^ 4S 


0.24 


T.89 




369- 


-1048 


(^FHT /FRPanOnS 


1 1 r\ T"i 't" 1 
l-^VJiicLLi 


0.06 




2453545.4 


2005.48 


0.24 


T.89 






510-715 


Skinakas 


Rcig 


2.0 




2453546.0 


2005.48 


0.24 


T.89 




369- 


-1048 


CFHT/ESPaDOnS 


Donati 


0.06 




2453547.0 


2005.48 


0.25 


T.89 




369- 


-1048 


CFHT/ESPaDOnS 


Donati 


0.06 


-0.14 


2453548.0 


2005.48 


0.25 


T.89 




369- 


-1048 


CFHT/ESPaDOnS 


Donati 


0.06 




2453554.6 


2005.50 


0.26 


T.89 


401- 


-590 


653-674 


OHP/Elodie 


Srv/Rauw 


0.3 


+0.37 


2453559.5 


2005.52 


0.27 


T.90 




391- 


-999 


Loiano 


Negueruela 


1.5 


+0.48 


2453563.4 


2005.53 


0.28 


T.90 




468- 


-716 


Skinakas 


Reig 


2.0 


+0.51 


2453565.5 


2005.53 


0.28 


T.90 






510-716 


Skinakas 


Rcig 


2.0 


+0.63 


2453568.5 


2005.54 


0.29 


T.90 


401- 


-590 


653-674 


OHP/Elodie 


Srv/Rauw 


0.3 


+0.66 


2453585.5 


2005.59 


0.32 


T.91 




390- 


-682 


OHP/Elodie 


Srv/Rauw 


0.1 


+0.94 


2453594.5 


2005.61 


0.33 


T.92 


401- 


-590 


653-674 


OHP/Elodie 


Srv/Rauw 


0.3 


+1.22 



18 /. D. Howarth et al. 



JD 


Year 


Phase <i>a 


Phase 




A range 


Telescope/ 


Observer 


AA 








(Ha) 


(orb) 




(nm) 


Instrument 




(A) 


(A) 


2453600 3 


2005 63 


0.34 


1.92 




362-715 






2.0 


+1.29 


2453601 5 


2005 63 


0.35 


1.92 




300—682 


OHP /Rlnrlip 


rtmiTP"!" /iNlPcmp'mpl a 


0.3 


+1.39 




9110^ R4 


0.35 


1.92 




R4R_RR9 


DAO 


i-JUllldlUCl 


0.3 


+1.50 




200^ 64 


0.36 


1.93 






DAO 


Rr^ n Ipn Hpt 

i— J WllldiLiCl. 


0.3 


-1-1 "^7 






0.36 


T.93 




R4S_Re4 


Tartu 




1.5 


+1.38 




900^ 6^ 


0.36 


1.93 




R47_RR9 


DAO 


n 1 oil Hpt 


0.3 


+1.30 


^'rtUOU J- tj .rt 


900"^ R7 


0.37 


T.93 




RQ7_R7Q 


Tartu 




1.5 


+1.38 


2453616.4 


2005 67 


0.37 


1.93 






Tartu 


Kolka 


0.5 


+1.54 


2453620.4 


2005 68 


0.38 


T.94 




640-662 


Tartu 


Kolka 


0.5 


+1.45 




XUUO.UO 


u.oo 


T 94 




RQ7_R7Q 


±al tU 


XVUlKcl 


-L.O 


+1 42 


Zl '± !j t>U Z • ^ 


9nn^ 70 


0.39 


1.94 


401- 


-5Q0 RR3— R74 


OHP /Plnrlip 


01 V / rXiCLLiW 


0.3 


+1.52 




9(10"^ 70 


u.oy 


T 94 




RQ7_R74 


±ai IjU 


iivoiita 


-L.O 


-Ul 4Q 




900"^ 71 


0.40 


T.94 




R4Q_RR9 


Tartu 


Knlka 

J:\.Ull\.cl 


0.5 


+1.38 




901^^ 71 


0.40 


T.94 




g37_g74 


Tartu 


KnlWa 


1.5 


+1.46 




900 77 


0.44 


T.96 


401- 


-'^QO Rfi'^— fi74 




01 V / lv,dU.W 


0.3 


+1.48 




900 77 


44 


T QR 


401- 


f^QO R"^*^ fi74 
oyu uoo u / 4: 


\J Lis: 1 JZjlOQKr 


S r V / Rau w 




U.O 


-Ul 

-p J- .00 




900"^ 77 


0.44 


T.96 




R4Q_RR9 


Tartu 


Knlka 


0.5 


+1.36 




900"^ CI 


47 


T Q7 




QOO— R>i9 


w n 1 1 cjioQie 


S r V / R,auw 


0*^ 

U.Uo 


-|- l.OO 


94^1fiR1 ^ 


900^ 


0.49 


1.97 


401- 


-5Q0 Ri^S— R74 
oyu uoo u 1 ^ 


OHP /Rlnrlip 


01 V / rXjCLLiW 


0.3 


+1.58 


94c;q7i q r 

Z'lOO 1 J-O.U 


900^1 


^R 


00 
u.uu 




R47— RR9 


DAO 


B ohleiicier 




U.O 


-Ul '^7 


94'=1'^799 fi 


900^ QR 


0.57 


0.00 




R47_RRO 


DAO 


J_> U 1 Udi 1 1 LlCi 1 


0.3 


+1.42 


94'=1'^777 1 


900R 1 1 


0.67 


0.04 




R4Q_R7Q 

u^o 010 


DAO 




0.3 


+1.03 


94'=1'^7J^9 1 


900R 1 9 


0.68 


0.04 




R4Q_RR9 


DAO 


UUliidiLlCl 


0.3 


+0.96 




900R 1 R 


0.71 


0.05 




R4q_RQ4 


Tartu 


KnlWa 


1.5 


+1.00 




900R 90 


0.74 


0.06 




R4Q_RQ4 


Tartu 


Knlkfl 


1.5 


+0.44 


94'^'^^^1 7 


900R 99 


0.75 


0.06 




R'i7_R7'^ 

uo / u / 


Tartu 


Knlka 

iVUllvcL 


1.5 


+0.14 


94'^'^S'^fi 
z^oooou .0 


900R 97 


0.78 


0.08 




'^S'=i— 707 
000 ( U i 


WHT /mm 


j-jdJiiiiuii 


0.6 


—0.68 


'± !j ^0 • !j 


9006 SO 


0.80 


0.08 




uo 1 U 1 


Tcirtu 


Kolka 


1.5 


— 1.18 




900R S4 


0.83 


0.09 




RQ7_R7Q 
UO 1 U 1 


Tartu 


Knlka 


1.5 


—1.87 


2453863 4 


2006 85 


0.83 


0.09 




637-673 


Tartu 


Kolka 


1.5 


—2.08 




900R S5 


0.84 


0.09 




RQ7_R7Q 
UO 1 U 1 


Tartu 


Knlka 


1.5 


—2.12 


94'=13SR7 4 


900R "^R 


0.84 


0.10 




RQ7_f;7Q 

uo 1 U 1 


Tartu 


Knlka 


1.5 


—2.09 


94'^'^K71 


900R "^7 


0.85 


0.10 




g37_g73 


Tartu 


Knlka 


1.5 


—2.32 




2006 "^0 


0.86 


0.10 




uo 1 U ( 


Tartu 


Knlka 


1.5 


—2.47 


2453889.5 


2006.42 


0.88 


0.11 




637-673 


Tartu 


Kolka 


1.5 


-3.00 


2453892.7 


2006.43 


0.89 


0.11 




386-705 


WHT/ISIS 


Lennon 


0.4 


-3.11 


2453895.4 


2006.44 


0.89 


0.11 




637-673 


Tartu 


Kolka 


1.5 


-3.19 


2453896.0 


2006.44 


0.89 


0.11 




369-1048 


CFHT/ESPaDOnS 


Donati 


0.06 


-3.17 


2453900.4 


2006.45 


0.90 


0.12 




637-673 


Tartu 


Kolka 


1.5 


-3.15 


2453941.5 


2006.56 


0.98 


0.14 




384-706 


WHT/ISIS 


Leisy/Lennon 


0.4 


-4.10 


2453955.4 


2006.60 


1.00 


0.15 




637-673 


Tartu 


Kolka 


1.5 


-4.00 


2453956.5 


2006.60 


1.01 


0.15 




637-673 


Tartu 


Kolka 


1.5 


-3.96 


2453958.4 


2006.61 


1.01 


0.15 




637-673 


Tartu 


Kolka 


1.5 


-4.14 


2453966.6 


2006.63 


1.03 


0.16 




384-706 


WHT/ISIS 


Leisy /Lennon 


0.4 


-4.07 


2454002.5 


2006.73 


1.09 


0.18 




384-706 


WHT/ISIS 


Leisy/Lennon 


0.4 


-3.32 



